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Abstract 

Using the method of correlation dynamics we investigate the properties of 
a field-theory for fermions and scalar bosons coupled via a Yukawa interaction. 
Within this approach, which consists in an expansion of full equal-time Green 
functions into connected equal-time Green functions and a corresponding trun- 
cation of the hierarchy of equations of motion we carry out calculations up to 4th 
order in the connected Green functions and evaluate the effective potential of the 
theory in 1+1 dimensions on a torus. Comparing the different approximations 
we find a strong influence of the connected 4-point functions on the properties 
of the system. 
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1 Introduction 



Up to now the problem of strongly interacting fermion and boson field-theories is not 
solved in a satisfying manner and there is a need for nonperturbative approximation 
schemes. In the last decades many non-perturbative approaches were developed such 
as variational methods |1|, 0|, Dyson-Schwinger-calculations p, |^ or the coupled clus- 
ter expansion ^. Among the approaches that are based on techniques known from 
many-particle physics is the method of correlation dynamics, which was successfully 
applied to the nonrelativistic nuclear many-particle problem H, || and for the de- 



scription of interacting boson systems [10|. Within a formulation in terms of connected 
equal-time Green functions the problem of ground-state symmetry breaking was stud- 
ied in detail for $^-theory in 1 + 1 and 2-1-1 space-time dimensions [|1T], |T^ as well 
as the convergence properties of the truncation scheme applied In this article the 
correlation dynamical approach is extended from a pure boson theory to a field theory 
containing both bosonic and fermionic degrees of freedom, i.e. Yukawa-theory in 1 + 1 
space-time dimensions. 

The paper is organized as follows: In Section 2 we present the connected Green 
function approach for Yukawa-theory in 1 + 1 space-time dimensions and discuss the 
correlation dynamical hierarchy for connected equal-time Green functions in different 
truncation schemes. Furtheron, we study the divergence structure of the Yukawa-theory 
in 1 + 1 dimensions and renormalize as in a comparable Gaussian-Effective-Potential- 
(GEP)-approach In Section 3 we show the results obtained in our numerical 

calculations for different truncation schemes up to the level of connected 4-point func- 
tions. A summary is given in Section 4 while the Appendices contain the higher order 
cluster expansions and all correlation dynamical equations of motion up to fourth order. 



2 Correlation dynamics of Yukawa-theory 

2.1 The equations of motion for connected equal-time 
Green functions 

In the present work we investigate the Yukawa-theory in 1 + 1 space-time dimensions 
given by the Lagrangian, 

C = U^^d^^-]^ml^^ + ^>{i^^'d^-MB)'^-gB^'^^. (1) 
or the Hamiltonian density, 

n = in2 + i(V<I>)' + im|$2 + ^(-z7V + MB)^ + (7B^^$. (2) 

Here rriB and Mb denote the bare masses of the boson and the fermion fields, respec- 
tively, while qb represents the bare Yukawa-coupling. The fields themselves have to 
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be considered as bare fields where we have dropped the corresponding indices. They 
fulfill equal-time commutation relations for the boson fields 

mx,t),U{y,t)] = zS^''\x-y), (3) 

mx,t),^y,t)] = [U{x,t),Uiy,t)] = (4) 
and equal-time anticommutation relations for the fermion fields 

{^.(f , t), My, t)} = 7°. S^''\x -y), (5) 

{^a{x,t),^b{y,t)} - {^a{x,t),^h{y,t)} = 0, (6) 

where a and b denote the components of the spinor. Making use of these relations one 
obtains the time evolution of any operator O via the Heisenberg equation 

tdtO = [d,H]. (7) 

The equations of motion for the field operators read: 

= U{x), (8) 

dtU{x) = {Vl-ml)^x)-gBMS)MS) , (9) 

dt ^a{x) = aiaV,^i{x) + iMePia^iix) + igBPia^i{x)^x) , (10) 

dt = - auV^^iix) - iMePu'^iix) - igePbi'^iixMx) , 

where all operators are considered at the same time t which is suppressed in our nota- 
tion. In 1-1-1 dimensions the Dirac-matrices are given by the repesentation 

For any arbitrary product of these field operators the time evolution is obtained cor- 
respondingly. The boson field momentum 11 is considered explicitly in order to end 
up with equations of motion of first order in time. Due to the Yukawa-coupling term 
the time evolution of any field operator product, which consists at least of two fermion 
fields or one boson field momentum, is influenced by fleld operator products of the next 
higher order. 

By taking the expectation value of these operator identities we end up with an infi- 
nite hierarchy of equations of motion for equal-time Green functions which is equivalent 



3 



to the BBGKYQ-density- matrix-hierarchy of ordinary many-body physics [|T5l . For the 
lowest orders we find: 

(<l>(xi)$(f2)) = {U{x^)^X2)) + {^x^)U{x2)) , (12) 

dt { U{xiMx2) ) = (V^^ - ml) { ^xiMx2) ) + ( U{xi)U{x2) ) (13) 

- 9b {^iixi)^i{xi)^ix2)) , 



(V^, - O ( Hx,Mx2) ) + {Vl - ml) { n(x 0$ (X2) ) (14) 
- gn ( ^i{xi)^i{xi)U{x2) ) ~ 9b { ^z(f2)^'z(^2)n(fi) ) , 

aia^cc^ ( ^i{xi)<i/b{x2) ) - auV^, ( ^a{xi)-^i{x2) ) (15) 

+ iMBpla { ^l{x,)^bix2) ) - iMbPu { ^„(Xi)^K^2) ) 
+ i9Bpla {^lixi)^b{x2)^iXi)) 

- igePbi {^a{xi)^ 1^x2)^x2) ) . 

In order to solve this infinite set of equations a truncation scheme must be applied. 
The ansatz of correlation dynamics consists of two steps: i) all full equal-time Green 
functions are replaced by sums of products of connected equal-time Green functions via 
a cluster expansion, which leads to a - still infinite - hierarchy for the connected Green 
functions, ii) a truncation scheme is applied to this hierarchy by neglecting all connected 
equal-time Green functions of order n > N . This closed set of equations of motion for 
connected equal-time Green functions is denoted as the correlation dynamical N -point 
approximation. 

Before deriving this set for the Yukawa-theory we separate the classical part from 
the boson field, since here we are primarily interested in the effective potential of the 
theory. The boson field is written as 

$ = $o + $, (16) 

where $o is the classical contribution given by ($) = $o and $ the pure quantum 
mechanical part with vanishing expectation value ($) = 0. With this separation the 
effective potential of the field theory is equivalent to the minimal energy density in the 
subspace of fixed vacuum expectation value i-e. 

$0, (S|S) = 1. (17) 



Mb - ^B$o)^ - 9b^^^ , (18) 



dt{n{x,)Uix2)) = 

dt{^a{Xl)^b{x2)) = 



V,s{%) = , with = 

Using (|16D the modified Lagrangian reads 

£ = ^9^$ 9^$ - im|(<l + $o)2 + ^{tYd^, - 
^Bogoliubov-Born-Green-Kirkwood-Yvon 
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and the corresponding Hamiltonian density is given by 

n = lfl^ + liV^f + Imli^ + $o)' + ^i-i^V + Mb + ^b$o)^ + gs^"^^ ■ (19) 

Since we take $o to be fixed, i.e. 9t$o = 0, we have 11 = dt^ = dt^ = 11. As a 
consequence the Hamihonian density now contains a contribution which is generated 
by the interaction part and has the same structure as the fermion mass term. The 
modified equations of motion for the full equal-time Green functions are obtained by 
replacing Mb by Mb + gs^a and the boson field $ by its quantum part $. 

The cluster expansion is derived from the relation between the generating function- 
als of full and connected Green functions and by taking the well defined equal-time 
limit. For the bosonic case this is presented in detail in Ref. in context of $^-theory 
and thus will not be repeated here. We only show the expansion for the lowest order 
boson Green functions (with field operators Oi{xi) G {$(xi), n(xj)} ): 

(6i(fi)) = (Oi(fi))., (20) 
(Oi(xi) 62(^2) ) = (Oi(xi) 62(^2) )c + (Oi(xi)), ( 62(^2) )c, (21) 

(6i(fi) 62(^2) 03(5^3)) = (Oi(xi) 02(^2) 03(5^3) )c 

+ (01(^1)02(5^2) )c (03(5^3) )c + ( Ol(^l) 03(^3) )c( 02(^2) )c 

+ ( 02(^2) 03(^3) )c ( Oi(xi) )c + ( Oi(xi) )c ( 02(^2) )c { 03(5^3) )c . (22) 

In the fermionic sector the source terms appearing in the generating functionals are 
anticommuting Grassmann fields. Thus the cluster expansions differ from the boson 
case and read in the lowest orders: 

(f(fl)^(f2)) = (^(fl)^(x2))c, (23) 

(^'(fl)^(f2)^(:?3)*(^4)) = (^'(fl)^(x2)^(x3)^(x4))c 
+ ( ^{Xi) ^(^4) )c ( ^(^2) ^(^3) )c 

- ( ^(fi) ^(f3) )c ( ^^(^2) ^(^4) )c . (24) 

Due to the absence of fermionic n-point functions of odd order the full and the con- 
nected 2-point function are equivalent. Therefore the index -c for the connected Green 
functions can be dropped for as in case of the boson 1-point function. For the 

mixed 3-point functions we find: 

(f(fl)^(f2) 03(5^3)}) = ('l'(Xi)^(f2)03(x3)}), 

+ (^(fi)^(f2))c( 03(5^3)); (25) 
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the higher orders up to the 5-point-level are given in Appendix A. 

An insertion of the cluster expansions into the hierarchy of equations of motion for 
full equal-time Green functions leads to a corresponding hierarchy for connected equal- 
time Green functions. In the present work three different truncation prescriptions are 
applied to this coupled set of equations, i.e. including all connected equal-time Green 
functions up to the 2-point-, the 3-point- and the 4-point-lcvcl. These approximations 
are denoted as Fi+iCD(2)-, Yi^iCD{3)- and Fi+iCD(4)-approximation. As examples 
for the equations of motion for connected equal-time Green functions we present those 
of lowest order: 

dt{Hxi)Hx2))c = {fl{x,)^X2))c + {Hx,)Il{x2))c, (26) 

dt { fl(f i) ^X2) )c = t{Xi) ( <l>(f i) ^X2) )c + ( n(fl) tl{x2) )c (27) 

- gB{^l{Xl)^l{Xi)^X2))c, 

dt { n(f l) fl(f2) )c = tiXi) { $(fi) fl(f2) )c + tix2) { fl(f l) $(f2) )c (28) 

- Qb { ^{Xi) ^'(fi) fl(f2) )c - 9b{ ^{x2) ^'(^2) n(^l) )c , 

dt { *„(fi) ^b{x2) )c = i{x^)la { t^f l) M^2) )c - i{S2)w { ^a{xi) ^1(^2) )c (29) 

+ igsPla { ^l{xi) ^b{x2) Hxi) )e 
- igsPbl { ^a{xi) ■^1{X2) I>(f2) )c , 

where the connected 3-point functions are neglected in the limit Yi+iCD{2). In (27)- 
(29) the following operators have been introduced for the bosonic part 

t{x) := (V^ - ml) (30) 

and the fermionic part 

i{x)ab ■■= aabV^ + i(Ms + C/B$o)/9a6 • (31) 

With regard to a numerical treatment we expand all field operators in a basis of or- 
thonormal single-particle wavefunctions. Here we use a set of plane waves in a box of 
size L with periodic boundary conditions. This choice allows for an explicit exploita- 
tion of the translational invariance of the theory. The expansion of the field operators 



is given by: 


l'(f) = 


J2 "Pec Xai^) > 

a 


n(f) = 


a 


^a(x) = 


a 


^b{x) = 


a 



(32) 
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where the elements of the single-particle basis Xa obey the following orthonormality 
relations in z/ + 1 space-time dimensions: 



d'^X Xa(^) Xp{.x) = 6af3 ■ (33) 

The greek indices refer to the basis elements and the latin indices represent the spinor 



components in (32). One straightforwardly obtains the equations of motion for the 
expansion coefficients of the connected equal-time Green functions. Introducing the 
abbreviations 



(a|AiA2) 



tap 



d^X X*a{x) 



B 



Xl3{x) 



d'^x x*a{x) t{x) X/six) 



7ab 



''aP - I d^X X*a{x) [aab^x + i{MB + gB%) (3ab\ Xfs{x) 

d^x x*A^) Lbix) Xpix) , 
they read in lowest order: 

dt{<fa<fl3)c = {T^af I3)c + {f 07^13)0 , 

X 

-9bJ2 ("1-^1-^2) (^iAi ^iA2 V5/3)c , 
A1A2 

X 

X1X2 
Ai A2 

dt { Ipaa i^bp )c = Yl ( ^bp )c " t^X ( ^aa i^nX )c 

X 

+ igBPlaYl ("l-^l-^a) (^ZAi V'6/3V'A2 )c 
Ai A2 

- igBPbl J2 (/^l^l^s) ( tpaa i^lXi ^X2 )c , 
Ai A2 



(34) 



(35) 



(36) 
(37) 

(38) 



(39) 



(40) 



Due to their length the equations of motion for further Green functions are shifted 
to Appendix B. A diagrammatical representation particularly for the time evolution of 
a mixed 3-point function is given in Appendix C. 

By using a finite basis of plane waves two cutoffs - ultraviolet as well as infrared - 
are introduced. For a given boxsize L the difference between two neighboring momenta 
is Afc = 27r/L such that a basis of n states in 1 + 1 dimensions contains an ultraviolet 
momentum-cutoff of kcut = {n — l)Afc/2. 



2.2 The renormalization scheme for Yukawa-theory 

Before carrying out any calculations on the basis of the correlation dynamical equations 
of motion one first has to renormalize the Yukawa-theory. In 1 + 1 dimensions the 
theory contains two divergent IPI-diagrams, the tadpole diagram and the polarization 
diagram (cf. Fig. |l]), as can be seen directly by evaluating the superficial degree of 
divergence. 

The tadpole diagram - due to its structure with only one external boson line - 
requires a renormalization of the external boson source. Therefore we have to consider 
a Yukawa-Lagrangian extended by an additional boson source term Jb The renor- 
malization is then done by introducing a source counterterm which cancels exactly the 
divergences arising from the tadpole diagram. In contrast to this prescription the bare 
source Jb is fixed in the GEP approach by the requirement that the energy density 
has to have a extremum at $o = (cf. [T^). The resulting difference between the un- 



renormalized and the renormalized source in GEP calculation is fully equivalent to the 
counterterm of the diagrammatical source renormalization. One has to keep in mind 
that this renormalization procedure for the boson source does not affect the correlation 
dynamical equations of motion of the connected equal-time Green functions. This re- 
sults from the fact that the boson source term only appears in the equations of motion 
for the 1-point functions. Since the 1-point functions are taken to be fixed by means 
of the separation ([T6|) no adjustments in the integration procedure are necessary. 

The handling of the divergent polarization diagram, however, is more sophisticated. 
It turns out that the only possibility to get rid of the divergence consists in a multi- 
plicative renormalization of the coupling constant. In the continuum the renormalized 
coupling constant has to be defined as 

gl = gl Jo(M) , (41) 

where Jq is an - in 1+1 dimensions logarithmically divergent - integral of the type 
dk 



In{rh) = J ^ ''(^) , ^k{m) = Vfc2 + m2 . (42) 

Thus for a finite renormalized coupling constant the bare coupling has to be infinites- 
imally small. We adopt the renormalization scheme of Ref. fl^ . 
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It is worth mentioning that there is no other renormahzation prescription possible. 
A mass renormahzation, which is the naive approach to a diagram of that structure, is 
not viable due to the equivalence of the renormalized masses and the effective masses 
at $0 = in the GEP-formulation. As a consequence the renormalized and unrenor- 
malized masses are the same for the bosons as well as for the fermions, respectively. On 
the other hand, a renormahzation of the wavcfunctions does also not lead to a solution 
since the effective potential will become unbounded from below. 

The GEP-calculation in the continuum leads to the following approximation for the 
effective potential of the Yukawa-theory in 1 + 1 dimensions: 



In order to compare the GEP solution with a correlation dynamical calculation one has 
to modify this expression due to the evaluation of the theory in discretized momentum 
space. In this case we remain with the following expression for a discretized GEP: 



In comparison to the continuum GEP the integrals are replaced by the corresponding 
sums: 



We note that in the GEP as well as in any correlation dynamical calculation the vacuum 

expectation values with respect to the perturbative vacuum are subtracted. Finally, 
the multiplicative renormahzation procedure for the Yukawa coupling constant has to 
be modified in the same way: 




(43) 





(45) 



gl = gl So{M) . 



(46) 
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3 Numerical results 



3.1 2- and 3-point approximations 

We first investigate the structure of tlie effective potential of the Yukawa theory in 
1 + 1 dimensions in the lowest orders of correlation dynamics, i.e. in 2-point and 
3-point approximation. In these calculations the effective potential is generated via 
an adiabatic switching of the coupling constant qr in time. The coupling constant is 
increased linearly - starting from the unperturbed groundstate with gR = Q - accord- 
ing to gR/m = at, where t denotes the propagation time and a the "adiabaticity" 
parameter. In the limit a ^ the Gell-Mann and Low theorem guarantees that the 



groundstate of the correlated system for fixed coupling qr is reached [|16|. For a de- 
tailed discussion of the applicability of this method we refer to Ref. where this 
has been investigated explicitly for $'^-theory. In the present case we only show that 
the applicability is given for all couplings in the correlation dynamical 2-point approx- 
imation. For the calculations we use m = M = 5 MeV while the system of equations 
is integrated in a basis of 19 plane waves and a boxsize L = 100 fm. 

To show the adiabatic convergence we display the real part of one coefficient of the 
connected Green function {ipiip2)c in yi_|_iCD(2)-approximation evaluated for various 
parameters a from a = 10~^ c/fm to a = 2000"^ c/fm for a bosonic field expectation 
value $0 = (cf. Fig. |^). With decreasing a one finds a reduction of the oscillation 
amplitudes. The corresponding increase of the oscillation frequency is a consequence 
of the presentation as a function of the dimensionless coupling gR/m and therefore an 
artefact of the time rescaling. For a < 500~^ c/fm the calculations have sufficiently 
well converged versus an asymptotic curve. All values of the observed quantity agree 
very well with a corresponding discretized GEP-calculation as presented by the solid 
line in Fig. ^. This also holds for all other coefficients of {^iip2)c- Furtheron - except 
explicitly mentioned otherwise - we use a = 2000"^ c/fm which ensures convergence 
in all coupling regimes investigated. In 3-point and 4-point approximation, however, 
this is no longer possible above some critical coupling (see below). 

In Fig. ^ the effective potential of the Yukawa-theory in 1-t-l dimensions is presented 
in Yi+iCD{2)- and for comparison in discretized GEP-approximation - as a function of 
the dimensionless coupling constant gR/m and the boson field expectation value $o- In 
both approximations one recognizes the parabolic structure of the effective potential 
for small renormalized coupling constants. With increasing coupling this parabola 
bends over and the effective potential turns negative above some critical coupling. In 
contrast to the pure parabolic structure of the continuum GEP (even for large gR/m) 
we obtain a double-well shape in the discretized calculation. This behaviour is a direct 
consequence of the ultraviolet cutoff induced by the finite number of single-particle 
basis states and can be explained by a Taylor-expansion of the discretized effective 
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potential (|^) around $o = 0, 



-m'-gl + glM'S.,{M) 



(47) 



Due to the ultraviolet cutoff the bare coupling constant does not become infinitesimally 
small but assumes a finite value. Therefore, in addition to the terms in the continuum- 
limit, further contributions of higher order in $o appear. Since the coefficients of the 
most important terms, the S'-i-terms, are positive, this leads to an increase of the GEP 
for high values of |$o|- 

Another feature of the correlation dynamical approximation on the 2-point-level 
within a finite set of basis states is the fact that for certain parameter regions the 
yi+iCZ)(2)-approximation differs from the discretized GEP-calculation. However, in 
the continuum limit with an infinite number of basis states both approximations gener- 
ate the same results since the GEP represents the stationary limit of the field-theoretical 
Hartree-Fock-Bogoliubov-(HFB)-approximation, which is equivalent to the correlation 
dynamical 2-point approximation. The difference in the discretized theory appears in 
all areas of the parameter plane where the effective fermion mass of the GEP-calculation 
M = M + qb ^0 becomes negative and can be interpreted as a result of a level-crossing 
effect. Beyond this hyperbola in the gn/m-^Q-plane the fermion single-particle state of 
lowest energy in the noninteracting case, i.e. the zero momentum mode, is energetically 
preferred in comparison to the highest lying fermion sea-state. Usually the highest sea- 
state and the first particle-state in Dirac's-picture are separated by twice the effective 
mass of the particle. In our case, where the effective mass becomes less than zero, 
this leads to an interchange of these states. Since in the Y'i+iCD(2)-approximation no 
pair-production in the zero momentum mode is possible within the adiabatic process, it 
does not reproduce the true groundstate on the 2-point-level in the parameter-regions 
where the level-crossing appears. 

In the Fi+iCD(3)-approximation we find two new features in comparison to calcu- 
lations on the 2-point-level. The first one is the existence of coupling regimes for small 
boson background fields where the propagation of the system of equations of motion 
breaks down. Thus results for the effective potential in Fi+iCD(3)-approximation are 
only accessible up to coupling constants of QR/m ^ 0.8 for all $0- In contrast, the effec- 
tive potential in 2-point approximation can be evaluated for all values of the coupling 
and the external background field. One finds that the correlation dynamical propaga- 
tion breaks down as soon as the effective potential assumes a non-convex structure (cf. 
Fig. ^) for the set of parameters {gfi/m,^Q). This effect is known from the correlation 
dynamical investigations of $^-theory and was studied in Ref . |jl3| . The reason for the 
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breakdown lies in the fact that the cluster decomposition property of connected Green 
functions is not valid anymore in a two-phase configuration as encountered in the "non- 
convex" region of an effective potential, where the correct energy density has a upper 
bound given via a Maxwell-construction |TB|, M that can be reached dynamically as 



soon as the approximation scheme allows for tunneling. Due to the inadequacy of the 
cluster decomposition the connected Green functions are not restricted to small values, 
but can dominate the full Green functions. In the latter case the correlation dynamical 
truncation is no longer a valid approximation. It is worthwhile to point out that this 
unstable behaviour in the present case is not a weakness of the correlation dynamical 
approach since for those coupling strengths where the discretized effective potential 
becomes non-convex, the system is no longer bounded from below in the continuum. 

The second aspect of the Yi+iCD(3)-approximation is the reduction of the total 
energy density. While the global shape of the effective potential - up to the critical 
coupling constant - is not changed much one finds a significant lowering of the energy 
density especially for boson background fields with small absolute values (cf. Fig. ^). 
This effect shows its maximum influence not for $o = 0, but for such combinations 
of the coupling constant and the external background field where the Yi^iCD{2)- 
approximation separates from the discretized GEP, i.e. where the effective fermion 
mass becomes zero. 

The reason for this behaviour is that in the li+iCD(3)-approximation (in contrast 
to the yL+iCD(2)-approximation) pair-production in the fermionic zero-momentum 
mode can take place. Thus the energy density can be lowered significantly in a region 
where such a state is energetically preferred. In order to understand why pair produc- 
tion in the fermionic zero mode is not possible in the Yi+iCD(2)-approximation, we 
treat the problem in particle-number representation. We consider the mode expansions 
(within a discretized momentum space) for all field operators 



$(f) 



a„ + a, 



k'i + rn? 



(48) 



nff) 



E 



(49) 




^g. Uia ~\~ Via 



(50) 



(51) 



where each annihilation operator annihilates the perturbative vacuum. 

Here the abbreviation -g for the inverse momentum index is introduced according to 



Va 



(52) 
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Evaluating the complete equation of motion for the fermionic particle-number operator 
we find (using Qq, = + aj^) : 

tdtibibj = ^bML ^ (llaAiA^) \ 

A1A2 Y 2 iJa iJAi ^^A2 

X {( ^l^Ai«A2 ) ^^aUi-X, + { &l4i«A2 ) ^^a^^Xl } 

-gBML J2 (l|«AiA2) ^ 

A1A2 Y 2 i^Q, i^Ai ^^A2 

X { ( ^Xi^a«A2 ) ^*Ai«ia + ( ^^Ai^a^Aa ) ^iX^^ia } ■ (53) 

In ( ^31) only operator products with 3 creation/annihilation operators lead to a nonzero 
r.h.s. (even for the zero momentum mode), such that for all momenta a particle- 
antiparticle production is possible. In a 2-point calculation where the expectation 
values of the operator products are factorized, only the following terms survive, 

{•a •p a^) = {a^) = V2mL $0(5^^5 . (54) 

In (^4]) the symbols represent an arbitrary fermionic operator with the corresponding 
momentum index. In a translationally invariant system one ends up with 

dtibibj = -2gB<^o^Im{{bidi)} (55) 

which becomes zero for the zero momentum mode. Thus when propagating the system 
in time - even with a simultaneous variation of the coupling constant - the explicit 
appearance of the momentum ka in (|55D prevents a change of the occupation number 
in the zero momentum mode. 



3.2 The 4-point approximation 

The highest order correlation dynamical approximation investigated in the present work 
is the Yi+iCD(4)-approximation. For these calculations we modify the size of the box 
to 5 fm due to the enormous increase in the number of differential equations that have 
to be integrated. In this approximation all connected boson 3-point functions as well 
as all connected 4-point functions, that contain fermion degrees of freedom, must be 
propagated. (Due to the special structure of the Yukawa interaction the pure bosonic 
connected Green functions of order 4 disappear in this approximation scheme.) For 
practical reasons we are limited to 11 basis states. One has to choose a smaller boxsize 
in order to achieve a sufficiently good ultraviolet convergence (at least for relative 
small boson background fields and up to coupling constants of gR/m = 0.8). All other 
parameters, that have been used in the previous calculations, remain the same. 
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In Y'i_|_iCi5(4)-approximation we observe a novel effect for small boson background 
fields, i.e. the propagation becomes unstable at $o = already for "low" coupling 
constants gR/m < 0.4. In the limit a one no longer achieves adiabatic convergence 
for gR/m > 0.4 (cf. Fig. ^). These coupling strengths are much lower than those for the 
non-convex phase of the system. To analyse this behaviour in more detail we introduce 
the concept of an equal-time self energy. One way to proceed this is to define a G- 
matrix G^^^ for the connected equal-time Green functions via the following relation 
between the full and the unconnected parts of the mixed 3-point function: 

E G'XX2 ( V'iAi ^bp ){^X,) = gs L^'^ ll E ( "l^l^2 ) ( V^iAi ¥^A2 ) . (56) 
A1A2 A1A2 

Due to translational invariance only the terms with kx^ = contribute on the l.h.s. of 
(|56D and the equation for the G-matrix reads 



A "^0 A1A2 



where 
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has been introduced. Defining the self energy by 

= 4f E G'Sa (Vx) = % G% , (59) 
we can write the equations of motion for the fermionic 2-point function as 

idt ii^aai^bp) = « E { ^Xa ( ^IX ^b/S ) " ^Ja ( V'aa ^^IX ) } 

A 

+ E { St ( ^IX i^bp ) - t'^X { ^aa ^Pix)} . (60) 



The self energies contain all contributions to the time evolution that exceed the 2- 
point-level. Keeping this in mind one can also define the self energies directly (and 
omit the singularity at $0 = in (p7D), such that the interaction part in the equation 
of motion has the same structure as the free contribution. With regard to the observed 
behaviour in Y'i+iCD(4)-approximation the self energy of the expectation value of an 
operator product of two fermionic particle operators is of interest: 

X X 

+ (?bMlE (1I/5A1A2) ^ 



A1A2 J'^^lS^Xi^X-z 
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(7bML ^ (l|aAiA2) ^ 



A1A2 Y 2 i^Ai t^A2 

X { ( 4i^/3«A2 ) ^*Ai«m + ( dx.bf^a^.^ ) t^ai^ia } • (61) 
Defining the self energy directly by 

^,A = E , ^ 

(^"^a) A1A2 ^2f]^l]AiWA2 

} , (62) 



we remain with the following expression for the equation of motion 

^dt{b%) = - E [^A. +txa] {b%) + E [^~/5A +S^a] {bib,) , (63) 

A A 

in which the free propagation part is given by the energy in the noninteracting case, i.e. 
iaf3 = fla^a/s- With the help of these self energies one can construct the renormalized 
mean-field 

= {Ret)^0 (64) 
and the renormalized mean-field-hamiltonian 

= iaP + KT- (65) 



The eigenvalues of ( |B5D are the interesting quantities because they correspond to the 
energies of the basis states in the interacting case. Plotting the energies versus the di- 
mensionless coupling constant gn/m we obtain the following picture (Fig. |^): Whereas 
the single-particle energies do not change much for gR/m < 0.2, at gR/m ^ 0.25 a 
substantial lowering of the energy-levels - especially in the zero-momentum mode - 
appears. For all other momentum modes this behaviour is not so strong and shifted to 
higher couplings gR/m. This contrasts the situation in the yi+iC-D(3)-approximation. 
In the correlation dynamical 3-point approximation the decrease of the single-particle 
energies is much slower. Here we have to mention that for this behaviour in the 
Yi+iCD(4)-calculation no - in comparison to the Yi+iCD(3)-approximation - addi- 
tional terms appear in the self energy. Only the magnitude of the terms in the self 
energy ( |5^ is changed significantly by the inclusion of the 4-point functions. Consider- 
ing the occupation numbers (b^b^) one finds - as a consequence of this rapid decrease 
of the single-particle energies - a very strong population of the corresponding modes for 
gR/m ^ 0.5 in the yi_|_iCD(4)-approximation (cf. Fig. j^). For small couplings the zero 
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momentum mode is occupied only weakly while the occupation number has its max- 
imum in the first mode; all other modes follow corresponding to their single-particle 
energies. When reaching the critical coupling regime the single-particle energy for the 
zero-momentum mode drops substantially such that it can be populated to a high de- 
gree. Finally, for increasing gn/m the occupation number of this mode exceeds those 
of all other momentum modes and the zero mode becomes occupied preferentially (cf. 
Fig. Here we also find a different picture in the Yi+iCi5(3)-approximation, where 
the lowest momentum modes are populated only for very high coupling constants and 
with much less intensity (cf. Fig. |]). 

The same behaviour can be seen in the alternative representation in Fig. |l^ (up- 
per part), where the occupation numbers of the modes as well as the dimensionless 
coupling constant are displayed logarithmically. For small coupling constants the oc- 
cupation numbers of all momentum modes are almost the same in both approxima- 
tions. The linear increase indicates a scaling of the occupation numbers with approx- 
imately the square of the (unrenormalized) coupling constant. We recall that, with 
exception of the zero-momentum mode, the occupation numbers of all other modes in 
Yi+iCD(3)-approximation agree quite well with those of the 2-point calculations. One 
can recover this dependence on QR/m analytically from the GEP-solution in the small 
coupling limit. Furthermore, one finds a scaling of the occupation numbers with $q 
if the absolute value of $0 is small. Thus for low gR/m the occupation numbers of 
all momentum modes (except the zero mode) are mainly determined by the 2-point 
contributions. The particle number in the zeroth mode is only influenced - due to the 
missing property of pair-production in the 2-point approximation - by contributions 
of higher order connected Green functions. Consequently the zero mode (although 
preferred by energy) is less populated than all residual modes. Its relative occupation 
number is increased in the yi+iCZ)(3)-approximation for growing coupling gR/m. In 
Yi+iCD(4)-approximation we find approximately the same particle number (compared 
to the Y'i+iCZ^(3)-approximation) for small gR/m. For coupling constants gR/m ~ 0.2 
the influence of the 4-point contribution begins to grow and first leads to an increase of 
the occupation number of the zero mode, then of all other modes with growing gR/m. 

We have to point out that the observed effect persists in the ultraviolet limit and is 
not an effect of the regularization. This we conclude from the fact that the decrease of 
the single-particle energies and the subsequent increase of the occupation numbers ap- 
pears in all momentum modes. We speculate that this new effect is due to condensation 
of fermion-antifermion pairs that only appears on the 4-point-level. 

In order to demonstrate that the steep rise in the occupation number of the zero 
momentum mode is not an artefact of the lower boxsize L we show in Fig. |10| (lower 
part) the occupation number of the lowest momentum mode as a function of gR/m for 
three boxsizes (L = 5 fm, 20 fm, 100 fm) in 3-point approximation using the usual 
parameterset of the 4-point approximation. One recognizes that for all boxsizes L 
the occupation number of the zero momentum mode increases only smoothly with 
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(66) 



gn/m and no breakdown appears in the displayed coupling regime. Without explicit 
representation we note that the occupation numbers for all higher lying modes in the 
3-point approximation essentially show a straight line as in the upper part of Fig. |I0|, 
however, shifted as a function of L. Therefore the difference between the Y'i+iCD(4)- 
and the Y'i+iCD(3)-approximation is substantial and not just an effect of the different 
L chosen in both approximations before. 

The regime of the classical boson field $0; where this behaviour can be found, is 
very small. Even for |$o| ~ 1-0 higher quantum correlations are suppressed so that the 
effect is no longer present. To illustrate the dependence on $o we show in Fig. 11 the 
correlation strength of the mixed 4-point-function (^f\E'nn), i.e. the quantity 

(which plays an important role for the applicability of the correlation dynamical ap- 
proach [^) as a function of QR/m for various $0- While the correlation strength 
increases smoothly for small couplings, it grows rapidly for a - $o dependent - critical 
coupling. Sharply peaked minima are a result of a change in sign in (^). As discussed 
in detail in a high correlation strength is allways connected with an increasing 
instability of the propagation which in the present case leads to a breakdown. Here 
we find that the stability of the system depends strongly on the absolute value of $o; 
for $0 = 0.0 the effect is observed first and is followed by $o-values of ±0.1 and ±0.2, 
respectively. For boson background fields $o ~ 1 the system becomes unstable only 
for relatively strong couplings gR/m ^ 1.5, which is in the regime of critical coupling 
constants as observed in the 3-point approximation. Thus the discussed effect in the 4- 
point approximation is limited to small background fields where quantum fluctuations 
are most prominent. 



4 Summary 

In the present work we have investigated the groundstate properties of the Yukawa- 
theory in 1 ± 1 dimensions using the method of correlation dynamics. We derived 
the infinite coupled hierarchy of equations of motion for connected equal-time Green 
functions and introduced three different truncation schemes, i.e. the 2-point-, 3-point- 
and 4-point approximation. Furtheron, we discussed the necessary renormalization 
prescription for the Yukawa-theory in 1 ± 1 dimensions which - besides an implicit 
counterterm renormalization of the boson source - consists in a multiplicative renor- 
malization of the coupling constant. 

In our numerical calculations we found the following results: 

• The generation of a correlated eigenstate of the interacting theory is possible by 
adiabatically switching-on the coupling constant (in line with the Gell-Mann and 
Low theorem). 
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The effective potentials in Yi^iCD{2)- and discretized GEP-approximation differ 
in the regime of "negative effective fermion masses" in case of a finite basis set 
due to an infrared- and ultraviolet-cutoff. These differences vanish in the infrared 
as well as in the ultraviolet limit and are due to the noninclusion of zero mode 
pair production in Y'i_|_iCD(2)-approximation. 

In the yi+iCD(3)-approximation the propagation breaks down in the non-convex 
region of the effective potential where the continuum theory becomes unbounded 
from below; this is also known from the correlation dynamical investigations of 



$ -theory and can be attributed to tunneling |13 . 



• In the Y'i+iCD(4)-approximation a novel effect is observed, which has no analo- 
gon in the pure bosonic theories investigated so far. Here the propagation breaks 
down already for small couplings gR/m, which is no longer a consequence of the 
transition to the non-convex parameter regime, but due to strong correlations 
induced by higher order connected Green functions in the fermionic sector. The 
associated lowering of the single-particle energies leads to a strong population 
of the previously weakly occupied zero momentum mode and finally to an in- 
creasing instability of the system. We speculate that this new effect is due to a 
condensation of fermion-antifermion pairs, but cannot proove this explicitly due 
to the limited set of basis states that can be employed numerically. 

The last effect, that persists in the ultraviolet limit, shows that novel phenomena 
appear in a theory with fermionic degrees of freedom since such a behaviour was not 
observed in any pure boson theory so far. Furthermore, it demonstrates the importance 
of quantum correlations of higher order especially for low or vanishing background fields 
$0- It remains to be seen if such effects will persist also in gauge field theories such as 



QCD where the correlation dynamical approach has been formulated in [20|. 
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A Cluster expansion of the mixed sector 



+ {^{xi 

+ {^{Xi 

+ {^{xi 

+ {^ixi 
+ {^ixi 
+ {^ixi 

+ {^{Xi 



(*(fl)*(f2 

-{^{xi 
-{^{xi 



03(X3) 04(^4)) = (^(fl)*(f2)03(f3) 04(^4) )c 
^(.^2) 03(^3) )e(04(.X4)) 
^(^2) 04(^4) )e(03(:?3)) 
^(^2))c (03(5^3) 04(^4) )c 
*(^2))c(C'3(f3))(04(x4)) 



03(:?3) 04(^4) 05(5^5)) = ('S'(^i) ^(5^2) 03(^3) 04(^4) 05(^5) )c 
^(52)03(:?3) 04(^4) )c(05(:?5)) 
^(52)03(.X3) 05(^5) )c(04(:?4)) 
^(5^2)04(f4)05(x5)),(03(f3)) 



^(5^2)03(f3)). 
^(fs) 03(X3) ), 
^(fs) 04(X4) ), 
^(fs) 04(:?4) ). 
^(fs) 05(:?5) 
^(f2) 05(^5)). 



( 04(^4) 05(x5))c 
(04(:?4))(05(:?5)) 
( 03(^3) 05(:?5))c 

{ds{xs)){d,ix,)) 
{dsixs)d,ix,)), 

(03(^3)) (04(^4)) 

^(^2) )c ( 03(^3) 04(^4) d^ix^) ), 

■^{X2))c{0s{xs)d,{x,)),{d,{x,)) 
*(^2))c (03(^3) 05(^5) )c( 04(^4) ) 
*(^2) >c( 04(^4) 05(^5) >c( 03(^3)) 
* (^2) >c( 03(^3)} (04(^4)) (05(^5)) 



*(f3)*(^4) 05(^5)) = (*(fl)^(f2)*(^3)*(^4)05(^5))c 
^(f2)*(f3)*(f4))c( 05(^5)) 
*(f4) 05(^5) )c(*(f2)*(f3))c 
*(f4)>c(^(^2)*(x3) 05(^5) )c 
*(f4))c(^(f2)*(f3))c( 05(^5)) 
*(^3) C'5(^5))c(^(^2)*(x4)>c 
^(f3))c(^'(f2)*(f4)05(f5))c 
*(f3))c(*(^2)*(f4))c( 05(^5)) 
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B Hierarchy of equations of motion for connected 
equal-time Green functions in a single-particle 
basis 



dt'\<^a<^ti^l)c = {T^a^l3^'y)c+ {^aT^fS ^-/) c + {VaVf3Tr-/)c (70) 

X 

- 9B^ {a\\l\2) {i^iX^i^iX^ VP V'y)c (71) 

A1A2 

C^t (TTaTT^ 9^7)0 = ( TT^ TT^ ) c + X! ( V^A TT/^ )c + ^/3A {T^aVxV-t)c 

X 

A1A2 

-9bY^ A1A2) ( V^iAi ^iX2 T^a )c (72) 
A1A2 

9t (TTaTT/? 777)0 = ^ taA {Vx'^0'^i)c + ^/3A {'^aVx'^j)c + tjX {'^a'^XV'r)c 
X 

- 9bY1 («I A1A2) ( IpiXi V'iAa % )c 

Ai A2 

- fl'B A1A2) ( i>iXi AX2 T^a 71-7 )c 

A1A2 

- 9bY1 (tI A1A2) ( Axi AX2 TTa TT^ )c (73) 

A1A2 

X 

+ WBPlaJ2 («|AlA2) {-4)1X^11)1,13 <PX2Vl)c 
A1A2 

- WbI^U X] (/?| A1A2) ( V'aa V'iAi </?A2 '/'t )c 

A1A2 

+ i9BPla J2 ("IA1A2) {i'lXii>b0)c {VX2Vl)c 
A1A2 

- i9Bl3bl A1A2) ( Ipaa V'iAi )c ( VX2 Vl )c 

A1A2 

+ {'f>aa'ipbf3'n-j)c (74) 
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A 
A 

+ igBplaYl ('^l'^l'^2) {i>l\i^bp^X2'^-y)c 
Ai A2 

- ^fi's/^W Yl iPl^i^^) { ^aa 'tplXi ^X2 T^-y )c 

A1A2 

+ i9BPla Yl («l^l^2) {'iI^ImAp)c (</?A2%)c 
Ai A2 

-igsPbl Y {lljaai>l\i)c {^X2^t)c 

A1A2 

-9bY (7I A1A2) ( ^aa ^iAi ^iA2 )c 
A1A2 

- 9bY (tI A1A2) ( -i/Saa V'iA2 )c ( 'i/'iAi -06/3 )c (75) 

A1A2 

dt{'l/jaa'i/jb0^y^s)c = Y ^Xa {i'lX^bp ) c " i^f^x{i>aa i>lX V^j ) c 

X 

+ igBplaY («l^l^2) {'ipl\i'ipbf3<PX2V'yVs)c 

A1A2 

+ i9Bpla Y {i>l\ii>bf}V^j)c {^X2^5)c 

A1A2 

+ igBPla Y ("l^l^2) {'4'l\^^b0^s)c {^X2'^-y)c 
A1A2 

- iQBpbl Y iPl^l^^) { ^aa IplXi VX2 ^1^5)0 

A1A2 

- WbI^U Y (/^I'^1'^2) ( ^aa V'ZAi 9^7 )c ( <^A2 <^5 )c 

A1A2 

-WbPuY {'4^aai^lXiV5)c {VX2Vi)c 

A1A2 

+ ('0aa'06/3 7r^</'5)c+ {i^aa^Pb^ ^-r '^s)c (76) 

dt{i'aa1pb0T^^^5)c = X) ^Aa {i'lXipbp ^s) c " ^Ja {i'aa'lplXT^^ ^s) c 

A 

+ ( V'aa TTa TT^ )c + X] ^tA ( V'aa V^A )c 

A 

+ igBplaY ('^l'^l'^2) {^l\ii>bf}VX2'^-yVs)c 
A1A2 

+ igBpla Y («l'^l'^2) {^l\ii'bf3T^'r)c {VX2^s)c 
A1A2 
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A1A2 

X1X2 

-igBpbl {lpaa'il^lXi7^j)c {<PX2Vs)c 

X1X2 

-igsPbl (/^l-^l-^s) {'ipaa'iplXi^s)c {^X2^'r)c 

X1X2 

-9bY1 (71^1^2) {i'aa1piXi^s)c {ipb/S 'lpiX2) c 
A1A2 

-9bY1 (7l''^l'^2) ( i^aa ^iXi )c ( V^fe/? ^1X2 T^5 )c (77) 
Ai A2 

dti'tpaalpblS'^j 7^6)0 = J2 *Aa {i^lxi^bH T^s) c " ^\ {i^aai^lX '^^'^5)0 

A 

+ Y ^7A {i^aa1pbpVx7^s)c + tsx {ipaa Ipb/d ^Px) c 
X 

+ igBPlaYl («l^l^2) {i^lXii^bpVX2'^lT^5)c 

X1X2 

+ igBpla Y {o^\^l^2) ('0ZAi'06/3 7r^)c ( V'Aa TTs )c 

Ai A2 

+ igBpla J2 ("1-^1-^2) {i>lXii^bp'^s)c {VX2T^l)c 
A1A2 

- igBl^bl Y (/5|^1^2) ( i^aa V'iAi VX2 ^ T^S )c 

X1X2 

-igBpbl Yl (/^l^l^2) {i'aa1plX^T^'^)c {^X2^s)c 
X1X2 

-igsPbl Y (/'^l^l^2) { i'aa'iplXi7^s)c ( V^Aa TT^ )c 
Ai A2 

- gsY (7IA1A2) {llJaaAxin)c {ApAx2)c 

X1X2 

-gsY (71-^1-^2) {i>aaAxi)c {i>b(3i>iX2T^s)c 
X1X2 

- gsY ^1^2) ( i^aa Axi % )c ( AX2 )c 

Ai A2 

- gsY (^l'^l'^2) ( V'aa V'iAi )c ( -06/3 V'iAa TT^ )c (78) 

Ai A2 

dt { i^aa i^bH tpds)c ^ Y ^Xa ( V'i A 'Ipb/S 'tpds)c + ( V'aa V'feA V'ct ) c 

A 

~ Y ^X {^aai'b(}Ax^d5)c + ^ii {^aai'bd ^c^ ^^ix) c 
X 
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C7 ¥'A2 /c 

A1A2 

A1A2 

+ igBPlbYl {'4>lM'^c-f)c {i'aa'il^d5<P\2)c 
A1A2 

-igsPwYl (/^l-^l-^s) {'ipl\i'ipds)c {i'aai'c-,^\2)c 
A1A2 

+ igBPcl (71^1^2) {i>aa1plXi)c {i^bfii^dS c 
A1A2 

-igBpcl (7IA1A2) {^bf3'iplXi)c {i^aa^dSVX2)c 
A1A2 

+ WBl3dl Y ('^l'^l'^2) {i'b0'4'lXi)c {i'aa^n^X2)c 
A1A2 

- «5'B/3(iJ X] ('^l^l^2) ( i'aa V'iAi )c ( V'b^ V'c7 </'A2 )c 
A1A2 
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C Diagrammatical representation of the equation 
of motion for a connected equal-time 3-point func- 
tion 

To clearify the structure of the correlation dynamical equations of motion we intro- 
duce a diagrammatical representation. In Fig. |12| we demonstrate this for the equation 
of motion of the connected equal-time 3-point function ( ipaa V'b/s % )c, where all contri- 
butions are displayed which arise from the ^- or the Il-part of the Green function (the 
contributions from field operator have the same structure as the ones from ^ only 
modified by a relative minus-sign). The n-point connected equal-time Greenf unctions 
are represented by triangles or rectangles with a corresponding number of field opera- 
tors while the greek letters at the lines - fermionic lines: solid, bosonic lines: dashed - 
going out of the connected Greenfunctions denote the momentum indices with respect 
to the chosen single-particle basis - spinor indices are neglected for transparency. We 
mention that all diagrams that emerge due to the application of the cluster expansion on 
the equation of motion for the corresponding unconnected 3-point function contribute 
which are connected in the usual (Feynman-) sense, while all disconnected diagrams 
cancel out when going over from the hierarchy for unconnected Green functions to the 
hierarchy for connected Green functions. 
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I 



Figure 1: The divergent diagrams of Yukawa-theory in 1 + 1 dimensions: the tadpole- 
diagram with superficial degree of divergence D = 1 (left) and the polarization-diagram 
with D = (right). 
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Figure 2: One coefficient of the fermionic 2-point function {ipiip2)c as a function of the 
dimensionless renormaUzcd couphng constant gR/m for various adiabaticity parameters 
a. For a — > one observes convergence of the Yi_)_iCD(2)-results towards the GEP- 
solution (sohd hne). 
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GEP 19 momentum states Y^^^CD(2) 

M = m = 5.0 MeV 
L = 100 fm 



Figure 3: Effective potential in yi+iCZ^(2)-approximation as a function of QR/m and 
$0 (right part) and tlie effective potential in the discretized GEP-approximation (left 
part). 
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Figure 4: 3-point contribution to the total energy density for various values of the 

boson background field $o in l^i+iCD(3)-calculation; for $o ^ —1-0 (upper part) and 
for $0 ^ —1.0 (lower part). The propagation breaks down for $o in the non-convex 
region of the effective potential. As "breakdown point" we specify those couplings at 
which either the propagation collapses or the energy-density of the 3-point coupling 
part starts to decrease rapidly. 
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Figure 5: Difference of the effective potential in Yi^iCD{3)- and in Yi+iCD(2)- 
approximation. 
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Figure 6: 3-point contribution to the energy density in Y'i+iC-D(4)-approximation 
for $0 = as a function of the dimensionless couphng constant gR/m for various 
adiabaticity parameters a. 
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Figure 7: Single-particle energies for the individual fermionic momentum modes in 
Yi+iC£)(3)- and Yi+iC£)(4)-approximation as a function of the dimensionless coupling 
constant gR/m. Note the massive decrease in the zeroth mode in the Yi+iCiP(4)- 
calculation at gR/m w 0.25. 
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Figure 8: The occupation numbers {bl^b^) of the zeroth and first fermionic parti- 
cle mode in Yi+iC£)(4)-approximation as a function of QR/m. In the corresponding 
Yi+iCD(3)-calculation the increase is only moderate and the occupation number of the 
zero-momentum mode remains small even for strong couplings. 
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Figure 9: Fermionic particle occupation numbers of the individual momentum 

modes as a function of the dimensionless coupling constant QR/m = 0.04... 0.64 (0.04) 
and the single-particle energy For coupling constants of about ~ 0.25 a strong 
reduction of the single-particle energies - especially of the zero- momentum mode - 
is observed. 
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Figure 10: Occupation numbers of the individual fermionic particle modes in 
Yi+iCD{3)- and in Yi_|_iCiI'(4)-approximation in the small coupling regime in double 
logarithmic representation (upper part). Occupation number of the lowest fermionic 
particle mode in Yi_|_iCD(3)-approximation for various boxsizes (L = 5 fm, 20 fm, 100 
fm) (lower part). 
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Figure 11: Correlation strength |(^*nn)c / (^^^1111)1 as a function of the dimen- 
sionless couphng constant QR/m and the classical boson field $o- 
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Figure 12: 
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Figure 12 ; Diagrammatical representation of the equation of motion for the connected 
equal-time 3-point function ( ipaa ^7 )c- Connected equal-time n-point functions are 
represented by triangles or rectangles with a corresponding number of field-operators. 
The greek letters at the lines going out from the connected Green functions - fermionic 
lines: solid , bosonic lines: dashed - denote the momentum indices with respect to the 
chosen single particle basis (spinor indices are neglected for clearness) . At the "vertices" 
momentum conservation is obeyed according to (0). Displayed are all contributions to 
the time evolution of ( ipaa '4'b/3 T^-y )c (part a)) that arise from its ip- (left diagrams) and 
TT-terms (right diagrams) (cf. ([75D). Parts b), c) represent contributions stemming from 
the free part of the Hamiltonian; the t-symbols correspond to the operators defined in 
(|35D and (^); Parts d), e) display the coupling to connected Green functions of order 
n+1 introduced by the Yukawa interaction; Parts f), g) show the coupling to products 
of connected Green functions of lower order as a result of the cluster expansion. 
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